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5 ■ Abstract 



C^. The Feynman-Schwinger representation provides a convenient framework 

for the calculation of nonperturbative propagators. In this paper we first 



K^ , investigate an analytically solvable case, namely the scalar QED in 0+1 di- 

mension. With this toy model we illustrate how the formalism works. The 
analytic result for the self energy is compared with the perturbative result. 
Next, using a x^^P interaction, we discuss the regularization of various di- 
vergences encountered in this formalism. The ultraviolet divergence, which 
is common in standard perturbative field theory applications, is removed by 
using a Pauli-Villars regularization. We show that the divergence associated 
with large values of Feynman-Schwinger parameter s is spurious and it can 



be avoided by using an imaginary Feynman parameter is. 
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I. INTRODUCTION 

In nuclear physics one is often faced by problems that require nonperturbative methods. 
The best known example is the problem of bound states. Even if the underlying theory may 
have a small coupling constant (such as in QED), and therefore allows the use of perturbation 
theory in general, the treatment of bound states are inherently nonperturbative. The n-body 
bound state is defined by the pole of the interacting n-body propagator. A perturbative 
approximation of n-body propagator does not produce the bound state pole location. This 
can most easily be seen by the following example: 

= IJ^X + X^+X^+... (1.1) 

1 — X 

Any truncation of the righthandside (perturbation theory) will fail to produce the pole 
which is located at a; = 1. Therefore it is essential that reliable nonperturbative methods 
that take all orders of interaction into account are developed. For this reason, numerous 
nonperturbative methods have been developed and succesfuUy used in the literature. Some 
of the best known examples are relativistic bound state equations [|I]-|^ , and lattice gauge 
theory [|[. 

Relativistic bound state equations provide a practical and intuitive framework to analyze 
bound states. However this practicality comes with certain drawbacks. In particular, the 
bound state equations in general lack gauge invariance. The second problem is associated 
with the fact that a completely self consistent solution of bound state equations is very diffi- 
cult. A completely self consistent solution requires solving infinitely many coupled equations 
for all n-point functions of the theory. Since this is an impossible task, one is either forced 
to model various vertices and interaction kernels or specify them perturbatively. 

The second and more recent approach is known as lattice gauge theory (LGT). LGT 
is a Euclidean path integral based approach which relies on the discretization of space- 
time. An economical lattice simulation with a small lattice size of 5^ requires roughly 
4 (links) X 8 (real SU(3) parameters) x 5^ (space — timepoints) = 20000 integrations. In 



general with larger lattice sizes this cost goes as 32N^. The disadvantage of discretization is 
twofold. The first one is the excessive computational time required for lattice simulations. 
The second problem is the anomalies caused by the discretization, such as the fermion 
doubling problem [Q. 

In this paper we discuss yet another method known as the Feynman-Schwinger Repre- 
sentation (FSR) 1^-0. Similar to the LGT, the FSR approach is also based on Euclidean 
path integral formalism. The basic idea in the FSR approach is to integrate out all fields 
at the expense of introducing quantum mechanical path integrals over the trajectories of 
particles. Replacing the path integrals over fields with path integrals over trajectories has 
an enormous computational advantage. In the FSR approach, a calculation similar to the 
example given above requires only 4N integrations, where A^ is now the steps a particle takes 
between the initial and final states. In addition to this enormous savings in computational 
time, the FSR approach also employs a space time continuum and therefore does not suffer 
from problems such as fermion doubling and the continuum limit. 

An additional motivation in studying the FSR approach is to determine which subsets 
of diagrams give the dominant contribution to the n-body propagator. This is particularly 
important in determining what kind of approximations are reasonable within the context of 
bound state equations. Therefore the FSR is a very promising tool to do this. In this paper 
we report on results for the 2-point function. In studies of hadronic systems like vrA^ one 
usually models the self-energy contribution through the lowest order one-loop graph. Also 
this is used as a starting point for the improved action procedure proposed by Lepage 0. It 
is clearly useful to compare such a lowest order approximation with the full results obtained 
from a FSR calculation. We study here as a toy model the scalar QED (SQED) in space 
and 1 time dimension, which can be solved analytically. The intriguing issue we also address 
in this paper is the difficulty found in the Euclidean action functional for the case of a 0^- 
theory. In applying the FSR to the 4-point function in the case of a (/)^-theory confined to 
generalized ladders one encounters a difficulty that the Euclidean action diverges. It was 
conjectured and demonstrated in a simple example by Nieuwenhuis [§] that this problem 



arises due to the continuation to the Euchdean metric. This problem is examined in detail 
here and we in particular show that there exists a regularization method to remove this 
divergence. As a result a clear prescription is given for dealing in a proper way with the 
Euclidean action in this case. 

The organization of this paper is as follows: In the next section we start by discussing the 
case of SQED in 0+1 dimension. This is a simple case which can be worked out analytically. 
We consider the one-body and two-body propagators. We compare the one-body result 
for the dressed mass with the perturbation theory result. In the third section we consider 
the case of scalar interaction x^0 in 3+1 dimension. We consider the issue of Wick rota- 
tion in Feynman parameter s, and present the FSR result for the one-body dressed mass 
obtained by Monte-Carlo integration. The result is again compared by the perturbation 
theory prediction. 

II. SCALAR QED 

Massive scalar QED in 0+1 dimension is a simple interaction that enables one to obtain 
a fully analytical result for the dressed and bound state masses within the FSR approach. 
In this section we compare the self energy result obtained by perturbative methods with 
the full result obtained from the Feynman- Schwinger representation. The Minkowski metric 
expression for the scalar QED Lagrangian in Stueckelberg form is given by 

CsQED = -m\' - \f^ + ]^^J?A^ - \\{dAf + {d, - ieA,)r{d^ + ieA^)x, (2.1) 

where A represents the gauge field of mass /i, and x is the charged field of mass m. We 
employ the Feynman gauge (A = 1). The presence of a mass term for the exchange field 
breaks the gauge invariance. Here the mass term was introduced in order to avoid infrared 
singularities which are present in 0+1 dimension. For dimensions larger than n=2 the 
infrared singularity does not exist and therefore the limit /i -^ can be safely taken to restore 
the gauge invariance. Since we confine ourselves to 0+1 dimension, the antisymmetric 



tensor F^^ vanishes. Therefore, in Euchdean metric and in 0+1 dimension, the scalar QED 
Lagrangian can be written as 



'-SQED 



m\^ + {dxf + \ /x'A^ + i {dAf + e\^A^ - teA{x*dx - xdxl 



(2.2) 



In preparation for the path integration which will be performed below, it is more convenient 
to cast the Lagrangian into the following form 

1 



rE 



SQED — X 



m' -d^ - 2ieAd - led A + e'A 



2 a2 



X + -A{^^''-^')A. 



(2.3) 



In order to construct a gauge invariant Green's function G, we introduce a gauge link f/(x, y) 



ie / dz A{z) 



U{x, y) = exp 

The two-body Green's function for the transition from the initial state $j 
X*{x)U{x,x)x{x) to final state $/ = X*iy)U*iy,y)xiy) is given by 



(2.4) 



G{y, y\x,x) = N jVx* jVxjVA $}$, e"^^ 



(2.5) 



where 



S, 



^ ^ ^SQED- 



(2.6) 



Performing the path integrals over x ^-nd x* fields one finds 



G{y, y\x, x) = N VA {detS)-'U{x, x)U*{y, y) [S{x, y)S{x, y) + S{x, x)S{y, y)] 



X exp 



--ld'zA{z){^'-d')A{z) 



(2.7) 



where the interacting propagator is defined by 



S{x,y) =< y 



1 



X > . 



(2.8) 



m^ — &^ — 2ieA d — iedA + e^A^ 
As in lattice gauge theory calculations we use the quenched approximation, detS* ^ 1. In 
order to be able to carry out the remaining path integral over the exchange field A it is 
desirable to represent the interacting propagator in the form of an exponential. This can be 
achieved by using a Feynman representation for the interacting propagator. 
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S{x, y) =<y\ ds exp -s{m^ - d"^ - 2ieA d - led A + e^A^ + it) 

Jo I 

where the s integration is along the imaginary axis. Let us now define 



X >, 



(2.9) 



U{x,y,s) =< y\exp 



-s(-d^ - 2ieA d - iedA + eM 



2 a2^ 



X > 



(2.10) 



where U{x, y, s) satisfies 
d 



ds 



U{x, y, s) = {d^ + 2ieA d + iedA - e^A^)U{x, y, s 



f2.111 



This is equivalent to the Schroedinger equation for imaginary time t = is, with Hamiltonian 



H{p,z) = {p + ieA{z)Y. 



(2.12) 



The matrix element of the interacting propagator Eq. ( [^.9D can be written in terms of a 
quantum mechanical path integral. We know from quantum mechanics that 



/POO 
Vqexp i C{q{t),q{t)) 

The Lagrangian for Eq. ( |2.12|) is given by 



dt 



(2.13) 



L{z,z) = iezA{z) 



(2.14) 



Therefore, the quantum mechanical path integral representation for this propagator is given 
by 



/•oo /■ 

S{x,y) = —i ds Vzexp 



is{w? + ie)-i/A I dTZ^{T) +ie j dr zA{z{t)) 



(2.15) 



where the boundary conditions are given by 2;(0) = x, z{s) = y. This representation allows 
one to perform the remaining path integral over the exchange field A. The final result 
for the two-body propagator involves a quantum mechanical path integral that sums up 
contributions coming from all possible trajectories of particles 



poo /"OO r r 

G = -j ds j dsj{Vz),yj{V 
where S[Z] is given by 



^ ixy ^ 



-S[Z\ 



(2.16) 



z(x) 




FIG. 1. A sample trajectory of each particle along with various interactions are shown. 

S[Z] = -iK[z, s] - iK[z, s] + V[z, z]. (2.17) 



The free and the interaction contributions to S[Z] are given as 



1 /■! 

(m^ + ie)s dr 



-H\z\ 



dz^ij) dz^^r) 



As Jo dr dr 

V[z,z] = -\ ® dr z{t) (j) df z{f) A{z{t) — z{f), ^), 

A(x,^) ^^^ ^"^ 



(2.18) 
(2.19) 
(2.20) 



27rp2 + ^2 2/x ' 

where A{x) is the interaction kernel. K[z, s] represents the mass term and the kinetic term, 
and V[z, s] includes the self energy and the exchange interaction contributions (shown in 
Fig. |1|). The contour of integration C follows a counterclockwise trajectory x ^ y ^ y ^ 
X — > X as parameters r, and f are varied from to 1. The self energy and the exchange 
interaction contributions, which are embedded in expression |2.19|, have different signs. This 



follows from the fact that particles forming the two body bound state carry opposite charges. 
The bound state spectrum can be determined from the spectral decomposition of the two 
body Green's function 



G(r) = Yl CnC 



-mnT 



(2.21) 



n=0 



where T is defined as the average time between the initial and final states 

T = -{yi + y^-x^-Xi). 



(2.22) 



In the limit of large T, the ground state mass is given by 



„„=„„, _^,„|G(T)l.-f^^«'l^l'-"^' 



T^oo dT 



SVZe-s[z] 



(2.23) 



A. The one-body case 

In order to be able to compare with the perturbative result later, we concentrate on the 
one-body case. The one-body propagator is given by 



G(0,T) = j ds f{Vz)oTexp 



iK[z,s]-Vo[z] 



(2.24) 



The integral of the self interaction Eq. ( 2.19 ) can analytically be performed 

Vo[z] = — I drzir) f df z(f) 6-"^'^^^''^^^^, 
4/i JO Jo 



2^2 



1 



-mT 



yUT 



(2.25) 
(2.26) 



where the boundary conditions were chosen as z{0) = 0, and z{l) = T. Next, the path 
integral over z can be evaluated after a discretization in proper time. Since the only path 
dependence in the propagator is in the kinetic term, the path integral over z involves gaussian 
integrals which can be performed easily by using the following discretization 



{V)oT^{N/4nsf/'Ul-,' dz. 



(2.27) 



The s integral can also be evaluated by saddle point method giving 



G(0, T) = iVexp 



—rriT "^ 



T e' 



:i-e- 



^^lT^ 



(2.28) 



2/i2 2/x3' 

This is an exact result for large times T. The dressed mass can easily be obtained by 
taking the logarithmic derivative of this expression. Therefore, the one-body dressed mass 
for SQED in 0+1 dimension according to the FSR formalism is given by 



M = m + 



2/x2 



(2.29) 



Simplicity of the SQED in 0+1 dimension also allows one to get an analytical result for the 
two-body bound state mass. It can easily be seen that the two-body result for the total 
mass is given by 
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+ 
1 



FIG. 2. Self energy for the SQED. 



e^ \ / e^ \ e^ 



AW. = [<n + -,j + (rn + ^ j " ^ = 2'n, (2.30) 

where the first two terms are due to the one-body contributions and the last term is due to 
the exchange interaction. The exchange contribution to the mass (up to a missing factor of 
two), was also reported earlier in Ref. 0. The interesting feature of the result in Eq. (|2.30|) is 
the fact that the positive shift of one-body masses are exactly compensated by the negative 
binding energy created by the exchange interaction. Therefore the total bound state mass 
is exactly equal to the sum of bare masses. Therefore in this simple case vertex corrections 
do not contribute to the bound state mass. 

In order to be able to compare with the FSR prediction Eq. (|2.29|) we consider the 
perturbative treatment of self energy. 

B. The perturbative result 

In this section we consider the perturbative treatment of the self energy and compare 
the perturbative result with the FSR prediction Eq. ( p.29|) . The self energy (Figure |) in 
0+1 dimension is given by 



oo 



2 I dk {2p — ky 2 f^ dk 1 



Performing the Wick rotation we get the following Euclidean expression 



np) = -e'l 



°° dk 



-oo 27r 
Evaluating the integral we find 



(2p - ky 1 



(P + /i2)[(p_A;)2 + m2] A;2 + /i2 



(2.32) 
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S(p) = - 



o2 r 



(z/i - 2py 



{im — pY 



^[m? + (iyU — pY] m[/i^ + (zTTi + pY 



The dressed propagator corresponding to this self energy is 

—i r .^ , n1 —i 



^^d{p) 



—I 



+ 



m^ _l_ p^ ^72"^ + p^ 

—i 



-iTjEip) 



m^ + p^ 



/i 



+ 



(2.33) 



The coordinate space form of the dressed propagator is 



(2.34) 



Mt)=l ^ 



^ipt 



Ne 



-Mt 



27T \m'^ + p'^ + Tie{p) , 
where M is the dressed mass and A^ is a normahzation factor. The dependence of M 
couphng strength e can be obtained from the solution of the on-shell condition 



M= Jm'^ + T.E{iM), 



(2.35) 
on the 

(2.36) 



which must be real if the dressed mass is to be stable. Therefore, for SQED, the equation 
determining the dressed mass takes the following form 



2 



(/i - 2Mf 



(m - Mf 



1 
+ - 



(2.37) 



.lilrn^ - (/i - M)2] m[/i2 - (m + Mf] /ij 
This perturbative result can be compared with the exact result Eq. ( |2.29| ) found earlier. 
In figure § we present the case oi m = ^ = 1 GeV. For small values of coupling strength 
e^ the perturbative and the full results converge. From the figure we see, that although the 
higher loop contributions cannot entirely be neglected they are of limited size, suggesting 
that in this case the lowest one-loop contribution may be a reasonable approximation for not 
too strong couplings. This is consistent with the results from Ref. in the case of SQED in 
2+1 dimension. The perturbative result develops a complex mass beyond a critical coupling 



■"crit 



0.343 (GeV)^. At the critical point a 'collision' takes place with another real solution 



of Eq. ( |2.37| ), leading to two complex conjugated solutions with increasing e^. This happens 
at M = 1.49 GeV. This is an inadequacy of the perturbative approach. 

The occurrence of complex ghost poles in the propagator has also been found in Lee-like 



models [TyjTlIl and in vr — A^-interaction models [0 . Moreover, it is also interesting to note 
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0.2 0.3 

e'(GeV') 



FIG. 3. The function M{g'^) calculated by the FSR approach and the perturbative methods 
for values of m = /i = 1 GeV. According to the perturbative result there is a critical point at 



Jcrit 



0.343 (GeV)'^ beyond which the dressed mass becomes complex. 



that a similar critical behaviour was also observed within the context of one body Dyson- 



Schwinger equation in Ref. |T3[. In the Dyson-Schwinger-Bethe-Salpeter studies of hadron 
structure the lack of real and finite mass poles in the quark propagator is usually considered 
as an indication of confinement. On the other hand, the simple example of SQED study 
in 0+1 dimension shows that while the exact result for the dressed mass obtained from the 
FSR approach produces a real mass pole for all values of the couphng, the simple bubble 
sumation leads to complex mass poles for large coupling values. Therefore the connection 
between confinement and lack of real mass poles is far from clear. 

Having presented the study of SQED in 0+1 dimension, where analytical results are 
easily obtained and compared with the perturbative ones, we move on to the scalar x^0 
interaction in 3+1 dimension. 
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III. THE FEYNMAN-SCHWINGER FORMALISM FOR SCALAR FIELDS 

We consider the theory of charged scalar particles x of mass m interacting through the 
exchange of a neutral scalar particle of mass /i. For this case the analytical integration of 
path integrals are not possible and one needs computational tools. 

The Euclidean Lagrangian for this theory is given by 

C-E = X*K -d' + 9<P]X + I Hf^' - d'^. (3.1) 

The two body Green's function for the transition from the initial state $i = x*{^)x{^) to 
final state $/ = X*{y)x{y) is given by 

G{y,y\x,x) = N J Vx* J Vx J V(j)<^}<^ie-^^ . (3.2) 

Performing the path integrals over x and x* fields one finds 

G(y,y|x,a;) = iv|p</)(det5)-i[5(x,y)^(x,y) + ^(x,x)5(2/,y)]e-^/'^'^<^('''-^')^ (3.3) 

where the interacting propagator is defined by 

S{x,y)=<y\— — , I ^ > • (3-4) 

m^ — a"^ + gcp 

In order to be able to carry out the remaining path integral over the exchange field (p it is 
desirable to represent the interacting propagator in the form of an exponential. 

S{x, y)=<y\ /"" ds e-'(m'-9'+94>+ie) I 3, > _ (35) 



/o 

Here we want to comment on a subtlety about the Feynman representation. In earlier works 
PJ^ the following Feynman representation was used 

POO 

S{x,y)=<y\ rfse-^(™'-^'+5*)|a;>, (3.6) 

Jo 

The validity of this representation depends on the sign of the field which can be either 
positive or negative. If one accepts this representation, the problem manifests itself as an 
exponentially diverging s dependence after the path integral over (p is performed. In order 
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to circumvent the problem of the s singularity we use the Feynman representation given in 
Eq. (|3.5|). Let us define 



U{x, y, s) =< y |e-^(-^'+3^) \x>, (3.7) 

where U{x, y, s) satisfies 

^^U{x,y,s) = id'-g<P)U{x,y,s). (3.8) 

This is equivalent to the Schroedinger equation for imaginary time t = is, with Hamiltonian 

H{p,z)=p'-g<P{z). (3.9) 

The Lagrangian for the Hamiltonian given in Eq. ( p.9|) is 

L{z,z) = ^ + g<P{z). (3.10) 

Therefore, the interacting propagator can be expressed as 



/•oo p r PS PS 

S{x,y) = —i ds Vzexp is{rn?' + ie) —i/A / dr z^{t) + ig / dr (f){z{T)) 
Jo J I Jo Jo . 



(3.11) 



where the boundary conditions are given by z{0) = x, z{s) = y. The final result for the two- 
body propagator involves a quantum mechanical path integral that sums up contributions 
coming from all possible trajectories of particles. The only difference from the SQED case 
Eq. (|2.15|) is the replacement of ezA{z{T)) -^ gcj). Therefore, for the two body propagator 



one arrives at the same expresion as Eq. ( p.l6|) except now the new definition (compare with 
Eq. ( p.l9|) ) of the interaction term is 

V[z, z, s, s] = Vo[z, s] + 2V,2[z, z, s, s] + Vo[z, s], (3.12) 

where 

Vo[z,s] = ^-s^ C dr f'dT'A{z{r)~z{T'),fi), (3.13) 

2 Jo Jo 

_ _ g'^ _ f^ (-1 

Vi2\z,z,s,s] = —ss dr dr A{z{t) - z{t), fi). (3.14) 

2 Jo Jo 
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Here the Vo[2;,s] term represents the self energy contribution, while the Vi2[z,z,s,s] term 
represents the exchange interaction (Fig. |1]). The notable difference compared to the SQED 
case is that the interaction terms now depend on the s variable. The second difference from 
the SQED case is the fact that self energy and exchange interaction terms have the same 
signs. The interaction kernel A(x) is defined by 

A(a;,^) = / ,f^, f[\ = -^Ki(/i|x|). (3.15) 



The time of propagation, T, is defined as before in Eq. ( |2.22| ) 

In principle one can work with equation ( |2.23| ), using the interaction given in Eq. ( p. 121 ) 



to determine the ground state energy of the bound state. However this is in practice very 
costly. The problem is the oscillatory behavior of the integrand which forbids the use of 
Monte-Carlo techniques for integration. Therefore it is desirable to make a Wick rotation 
in variable s. In the next section we discuss how this rotation can be made without leading 
to a large s divergence associated with the interaction term. 



A. The large s behavior and Wick rotation 

The one body propagator is given by 
G = i ds Vz exp 
where the s- independent functionals k'^[z] and v[z] are defined by 



A-2 

2 • 2 

tm s — I s V — es 

As 



(3.16) 



.2r.i_ ['dz,{T)dz'^{T) 



e[z]^ / dr-^— ^, (3.17) 

JO dr dr 



v[z] ^^ f dr f dr'Aizir) - z(r'),/i). (3.18) 

2 JO JO 

The path integral is discretized using 

{V),y ^ {N/ArrsY'^Ul-l' J d'z,, (3.19) 

where the s-dependence is critical in obtaining the correct normalization. The one body 
propagator after this discretization is given by 
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iV \ / .._^ /\/- 1 , / as 2 2 

im s — I- s V 

As 



G^^Q /n-..f Jiexp 



(3.20) 



This is a well defined and finite integral. In principle the number of steps A^ should be 
taken to infinity. If one keeps the N fixed, a simple replacement oi s —^ is clearly leads to 
a divergent s integral and is therefore not allowed. In order to put this integral into a form 
that allows Wick rotation without changing the physics we use the following trick. At large 
values of s the integral in Eq. (|3.20|) is highly damped because of the v and s^^ terms. The 
integrand is also highly oscillatory as s — > 0, or s ^ oo and therefore these regions do not 
contribute to the integral. In the limit f?^ — > the dominant contribution to the integral in 
Eq. ( |3.20|) can be shown, by using the saddle point method, to come from 



k 
s = iso = i-—. (3.21) 

2m 

Since the large s values do not contribute to the integral even without the interaction term, 
it is a good approximation to suppress the g"^ term at large s values. While this suppression is 
done it is important that the integrand is not modified in the region of dominant contribution 
s ~ isQ. This can be achieved by scaling the s variable, in the interaction term only, by 

(3.22) 



R{s,so)' 
where 

R{s,so) = l-{s-isof/T\ (3.23) 

In the free case, {g"^ = 0), the width W of the region of dominant s contribution goes as 



Therefore, in the free case the dominant contribution to the s integral comes iTom.i{sQ—W) < 
s < i{so + W). In order to ensure that the scaling given in Eq. ( |3.22D does not make a 



significant change in the region of dominant contribution, F should be chosen such that 

T>W (3.25) 
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(r/W) 



FIG. 4. The figure shows the insensitivity of the dressed mass to the width T for T >W. This 
case represents g = 5 GeV, and mT = 40. Results were obtained averaging about 4 Monte-Carlo 
runs at each T. 

It should be pointed out that as one increases the coupling strength, the value of Sq will 
deviate from its free value. Therefore, in general, sq has to be defined self consistently by 
monitoring the peak of the s distribution. In Figure (^) we display the insensitivity of the 
dressed mass to the width F, ioi g = 5 GeV and mT = 40. The results we present in the 
remainder of the paper are obtained with a choice of F^ = 2W'^. As a result of the scaling 
Eq. ( ^.22| ) the interaction term disappears at large s values where the integrand does not 
contribute anyway. The benefit of this replacement is in the fact that even though N is kept 
finite one can perform a Wick rotation rigorously in variable s to find a nonoscillatory and 
finite integral. Now let us take a closer look at the Wick rotation. 

After the redefinition given in Eq. (|3.22|) the s integral in Eq. (|3.16| ) takes the following 
form 
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Im(s) 




(0,0) 



Re(s) 



FIG. 5. Wick rotation 



f^ ds 
(j oc / ^777 exp 



s 



2N 



2 • "" 2 2/ Ti2 / \ 

im s — I OS K [s] V 



(3.26) 



The exponent has a singularity in the complex s plane at Sp = isQ ± T. One of these 
singularities is on the path of the Wick rotation. However it can easily be seen that it 
does not contribute to the integral. In particular, the contribution of the singularity at 
Ss = iso + r, let's call it Gpoie, is given by 
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(3.27) 



which is identically equal to 0. The contribution to the integral coming from the quarter 
circle can also be shown to vanish. On the quarter circle the interaction term approaches 
zero, and the integral is dominated by 1/ s^^ term which vanishes as the radius of circle goes 
to infinity. Therefore the integral vanishes along the quarter circle. The vanishing of the 
interaction term on the quarter circle is only possible if one assumes that the radius of the 
quarter circle is much greater than F. Therefore F can not be taken to infinity until the s 
integral is performed. 

Thus, one can indeed perform the Wick rotation without any complication to find a finite 
and nonoscillatory expression for the fully interacting two-body propagator: 



G = j dsj ds jiVz)^yjiVi 



Ixy 



xexp 
where 



-K[z, S] - K[z, S] + Vo[z, Sr] + Vo[z, Sr] + 2Vi2[z, Z, Sr, Sr 



(3.28) 



The discretized versions of kinetic and interaction terms are given by 

N ^ 
K[z, s] ^ (m^ + te)s - — J2{z, - z^.^f, (3.30) 

^^ i=i 
g'^s^ ^ 1 

M^^s] -^ 1^ Y. ^iij^i^i + ^i-i - % - Zj-l),^l), (3.31) 

Vi2[z,z,s,s\^ -— J2 A{-{zi + Zi_i-Zj-Zj_i),fi). (3.32) 

Having outhned the treatment of the large s behavior and the Wick rotation, we next address 
the regularization of the uhraviolet (short distance) singularities. 



B. The ultraviolet regularization 

The ultraviolet singularity in the kernel A(x, /i) Eq. ( p.l5|) can be regularized using a 
Pauli-Villars regularization prescription. In order to do this one replaces the kernel 

A(a;,/i) — > A{x, fi) — A{x,afi), (3.33) 

where a is in principle a large constant. The ultraviolet singularity in the interaction is of 
the type 

f dzzA{z,fi). (3.34) 

At short distances the kernel A{z,fi) goes as l/z"^. Therefore, we have a logarithmic type 
singularity. The Pauli-Villars regularization takes care of this singularity. The Pauli-Villars 
regularization is particularly convenient for Monte-Carlo simulations since it only involves 
a modification of the kernel. In order to achieve an efficient convergence in numerical 
simulations we use a rather small cut-off parameter a = 3. This choice leads to a less 
singular kernel. However this is not a major defect since the value of a can be increased 
arbitrarily at the cost additional computational time. 
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FIG. 6. The self energy for x^*/* interaction. 
C. Perturbation theory result for self energy 

In this section we study the self energy of a particle of mass m in lowest order in pertur- 
bation theory. We carry out the study in 3+1 dimension. 

The lowest order "bubble" diagram is shown in Fig. p. In 3+1 dimensions this diagram 



IS 



Eip)=tg'l 



d^k 



1 



(3.35) 



(27r)4 (m2 -{p- kf) (/i2 - P) ■ 
In order to compare the perturbative result with the FSR prediction we use the same 
regularization method, namely the Pauli-Villars regularization. Assuming that p^ < {m + 
/i)^, the integration over k may be rotated from the real axis to the imaginary axis without 
meeting any singularities. Substituting k -^ —ik and p —* —ip and gives the Euclidean 
expression for the self energy 

d'k A^ - ij^ 



^e{p) = -9'J 



(3.36) 



(277)4 (^2 + (p _ ^)2) (^2 + ^2) (^^2 + p) ' 

where the Pauli-Villars regularization mass A is chosen to be A = 3/i. Using the Feynman 
trick the integral can be evaluated giving 



Se(p) 
where /(p^,/i) is defined by 
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/(p^/x)-J(p^A) 



(3.37) 



J(p^,/i)= / (ixln[m^x + p^x(l — x) + /i^(l — x)], 
Jo 
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(3.38) 
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where 

The dependence of M on the couphng strength g can be obtained by analytic continuation 
of the Euchdean form of S given in Eq. (p.37|) . It is found from the on-shell condition 
Eq. (|;36D that 

,2 



ji^2 2 9 

M = m 



167r2 



/(-M',/i)-/(-M',A) , (3.40) 



The mass M is therefore the solution of this equation which must be real if the dressed mass 
is to be stable. 

The FSR result is obtained through Monte-Carlo integration. The dressed mass nio = 
liniT^^oo M{T), which is given by Eq. (|2.23|) becomes largely independent of T at large times 
mT > 40. As the coupling strength is increased the plateau is shifted towards higher T 
values. In Figure ^ we demonstrate how the stability is achieved as T increases for the case 
of ^ = 5 GeV. 

The behavior of dressed mass M{g^) as a function of the coupling constant is illustrated 
in Fig. ^. M{g'^) is always smaller than unity, and decreases as g increases. The agreement 
of the FSR result with the perturbative prediction is very good at low g"^ < 30 (GeV)^. We 
see that the mass shift is negative, corresponding to an attractive interaction. This should 
be contrasted with the SQED, where a positive mass shift is predicted. From the figures we 
see that the higher loop contributions increase the mass shifts in both cases. 

Moreover, the perturbative result displays a critical point as in the 0+1 dimension SQED 
case. According to the perturbative result (see Fig. |^), the dressed mass M decreases up to 
a critical critical value gcrit which occurs when the mass reaches to M^-it = 0.094 GeV. For 
this simple case the critical coupling is given by g^rit = 22.2 GeV. For larger values oi g there 
are no real solutions, showing that the dressed particle is unstable. For g > gcru the state 
does not propagate as a free particle. For the example shown in the figure, m = 1 GeV, and 
/i = 0.15 GeV. 
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FIG. 7. The Monte-Carlo result for the function M(T) is shown. Error bars reflect statistical 
errors associated with the Monte-Carlo sampling. The plateau region occurs around T > 40. 
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FIG. 8. The FSR result for the function M{g^) obtained from the Monte-Carlo simulation is 
shown along with the perturbative result. Error bars on the Monte-Carlo result are due to the 
fluctuations of the correlation function (which goes as e""**) under time. 
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FIG. 9. The perturbative result for the function M{g^) is shown for m = 1 GeV, and fi = 0.15. 
The critical coupling is given by ^crit = 22.2 GeV. For larger values of g there are no real solutions 

IV. CONCLUSIONS 

In this paper we have considered the SQED interaction in 0+1 dimension and the scalar 
X^<f) interaction in 3+1 dimension. We have shown that for the SQED, the analytical FSR 
result and the perturbative one are in agreement at small couplings. The exact SQED result 
for the dressed mass is real while a lowest order bubble sum produces complex mass poles as 
the coupling constant is increased. This example shows that the lack of real and finite mass 
poles do not necessarily imply confinement unless they are obtained by fully nonperturbative 
calculations. 

For the x^0 interaction we have shown that it is possible to perform a Wick rotation in 
Feynman parameter s and obtain a convergent expression for the s integration. 
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